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We propose a new ground state trial wavefunction for a two-dimensional Wigner crystal in a strong 
perpendicular magnetic field. The wavefunction includes Laughlin-Jastrow correlations between 
electron pairs, and may be interpreted as a crystal state of composite fermions or composite bosons. 
Treating the power m of the Laughlin-Jastrow factor as a variational parameter, we use quantum 
Monte Carlo simulations to compute the energy of these new states. We find that our wavefunctions 
have lower energy than existing crystalline wavefunctions in the lowest Landau level. Our results 
are consistent with experimental observations of the filling factor at which the transition between 
the fractional quantum Hall liquid and the Wigner crystal occurs for electron systems. Exchange 
contributions to the wavefunctions are estimated quantitatively and shown to be negligible for 
sufficiently small filling factors. 



PACS numbers: 73.40.Hm, 73.20.Dx 
I. INTRODUCTION 

It was first argued by Wigner in 1934 that a system 
of interacting but otherwise structureless electrons can 
have crystalline jprder in the limit of low density and 
low temperatureo. The first experimental evidence of 
the Wigner crystal (WC) was found well over 40 years 
later in a two-dimensional i42D) system of electrons ad- 
sorbed on a helium surfaceo. Nowadays, semiconductor 
heterojunction devices are considered a very promising 
environment for observing the WC. The advantage of 
heterojunction structures comes from the fact that the 
2D electron plane is spatially separated from the donor 
layer, so that the influence of these impurities can be 
substantially reduced compared to bulk semiconductor 
environments. Furthermore, it is now well known that a 
strong magnetic field perpendicular to the 2D plane can 
effectively localize electron wavefunctions while keeping 
the kinetic energy controllecoO. Since this lessens the 
otherwise severe low-density condition, it is believed that 
the WC can be stabilized in a sufficiently strong magnetic 
field. 

On the other hand, the fractional quantum Hall (FQH) 
liquid is known to be the ground state in certain ranges 
of strong magnetic field&Q. In this strongly correlated 
liquid state, the Hall resistivity p xy is quantized at dis- 
crete values and the diagonal resistivity p xx vanishes at 
zero temperature. In contrast, p xx presumably diverges 
at T = in the pinned WC. The FQH effect for v = l/m 
(m odd) is now fairly-well understood in terms of the 
Laughlin wavefunctionB. 

In several recent experiments with high mobility sam- 
ples, a sharp phase transition from the FQH state to an 
insulating state was observed as the magnetic field was 
increased both in electroroEl and holeO systems. Some 
properties of the insulating state such as the activation 
gap in charge transport closely resembles those of pinned 
charge density waves, supporting the interpretation of 
this insulating state as a WCL, Theoretical calculatioas. 
of both the FQH liquid energyll3 and the WC energyM'til 



are also in good agreement with the experiments as to the 
critical value of the magnetic field at which the transition 
occurs for a given electron density. 

However, not all experimental findings of the insulat- 
ing state are consistent with the conventionally accepted 
theoretical understandings of the WC. (I) First, there is a 
discrepancy in the energy of charged excitations. Partic- 
ularly, transport experimentsH reveal that the activation 
gap is an order of magnitude smaller than the theoret- 
ically, estimated energy to create a point defect in the 
Wdi3ii3. (II) Moreover, even deep into the insulating 
phase, anomalous behavior is observed when the filling 
factor v is an inverse odd integer, which may be related, 
to the FQH effect. Specifically, transport experimentsE£l 
exhibit a dip in the diagonal resistivity p xx of the in- 
sulating state pear v = 1/7. Also, photoluminescence 
experimentai3~E3 exhibit structure near odd denomina- 
tor filling factors down to 1/11, which looks very similar 
to structure seen at higher fillings wkeije.|the FQH ef- 
fect occurs. (Ill) Finally, experimentsEJ'Eaca show that 
the Hall resistivity p xy in the insulating phase saturates 
at its classical value B/nec, just as in the FQH liquid 
phase. This behavior cannot be understood in terms of 
a model of thermally activated point defects that are es- 
sentially non-interactingEl Interestingly, (II) and (III) 
suggest that some characteristics of the FQH effect are 
shared by the insulating state. 

The unusual behavior of p xy has led to speculation 
that the insulating phase is not a WC at all, but rathe 
disorder-dominated state called the "Hall insulator" c3 
However, it has been shown that interstitial defects in a 
WC can also lead to Hall insulating behavior if one intro- 
duces Laughlin-Jastrow correlations between the intersti- 
tials and the lattice electronsOta. The correlation was 
found to lower the energiz-to create such defects. How- 
ever, more careful studieso of the above interstitial state 
using Monte Carlo simulations suggest that in order to 
obtain such a small excitation energy as found in experi- 
ment, one must introduce Laughlin-Jastrow correlations 
into the ground state as well. In this paper we explore 
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the energetics of ground state wavefunctions of this form. 

In what follows we will introduce trial wavefunctions 
that take the form of a Laughlin-Jastrow factor multi- 
plying a properly (anti-)symmctrizcd product of single 
particle states. The wavefunction introduced here, thus 
corresponds to composite fermion or boson statesEHllEl. 
The energies of these states are computed using quan- 
tum Monte Carlo simulations, and it will be demon- 
strated that such states are generically lower in energy 
than other lowest Landau level WC states in the liter- 
ature. Our computational method in its simplest form 
ignores exchange corrections; i.e., the state multiplying 
the Laughlin-Jastrow factor is approximated as a simple 
product of single particle states. An in- principle exact 
computational scheme in which permutations of the sin- 
gle particle states are sampled shows that this is an ex- 
cellent approximation, provided the filling factor is not 
too large. 

This paper is organized as follows: The trial many 
body wavefunction of the Laughlin-Jastrow-correlated 
WC is introduced and some of its properties are discussed 
in Sec. 0. In Sec. HI, the ground state energy is com- 
puted using a Monte Carlo simulation. Various aspects 
of the results are also discussed. Sec. IV is devoted to 
discussions on the effect of the exchange energy and va- 
lidity of our approximation. Finally, we summarize the 
findings in Sec. M. Some technical details of the energy 
calculation can be found in the Appendix. 



II. TRIAL WAVEFUNCTION 

The Hamiltonian of 2D electrons moving in a magnetic 
field B and interacting with the Coulomb energy is writ- 
ten 
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where and Pi are the 2D position and momentum of the 
i-th electron, and A is the vector potential from which 
the magnetic field is given by B = V x A. We will ignore 
the spin degree of freedom assuming that it is completely 
polarized by the strong magnetic field. 

Since our trial wavefunction is closely related to 
the Hartree-Fock wavefunction of the WC suggested in 
Ref. [32], it is worthwhile to summarize the aspects of the 
Hartree-Fock wavefunction here. Explicitly it can be 
written 



(2.2) 



where (f>^ F are single electron wavefunctions and A an- 
tisymmetrizes the total many body wavefunction. To a 
first approximation in the large magnetic field, all elec- 
trons will lie completely in the lowest Landau level. The 
un-normalized single particle wavefunction is thus given 
by 



This wave function describes an electron localized at 
within a 2D Gaussian wavepacket. The magnetic length 
^/hc/eB determ ines the size of the wavepacket. The 



phase factor in Eq. (2.3) ensures that <fi^[ F is a product 



of e~' r '/ 4 ^ and an analytic function of z = x + iy, mak- 
ing it entirely lie in the lowest Landau level. In the large 
B limit, £ — ► 0, and the electrons become highly local- 
ized, behaving almost like classical point charges. The 
kinetic energy is given by the lowest Landau level energy 
Tiuj c /2 = heB/2mc, and is the same regardless of the 
Gaussian center R^. This allows one to treat R^ as vari- 
ational parameters in minimizing the total energy with 
respect to the Coulomb interaction. For a classical sys- 
tem, a triangular lattice is well knows to have the lowest 
Coulomb energy for a given densityEj. Therefore, in the 
limit B — > 00, the ground state is expected to be rep- 
resented by the above wavefunction, with R^ forming a 
triangular lattice. 

In a finite magnetic field, quantum fluctuations around 
the lattice sites become important and the above classical 
analogy is only approximate. Consequently, ^hf is not 
guaranteed to give the lowest energy at finite B. How- 
ever, if B is large enough, ^fjjF is still very close to the 
true ground state. For this reason, it has been used even 
at finite B by many authors, producing very good results. 
In this paper, however, we will impro\p,upon it by intro- 
ducing a correlation. Previous studiealil have introduced 
correlation factors that are exact for a harmonic Hamil- 
tonian; however, such wavefunctions allow fluctuations 
in which particles may occasionally closely approach one 
another. A correlation factor which by now is well-known 
to suppress such fluctuations is the Laughlin-Jastrow fac- 
tor: ni<j( z i — z j) m ') where z% = x% + iyi is the complex 
notation of the electron coordinates. Because of the ex- 
tra phase accumulated when one particle encircles a sec- 
ond, wavefunctions of this form may be understood as 
being comprised of particles that have m magnetic flux 
quanta attached to them. The idea of constructing wave- 
functions of this general form was first suggested in the 
context of the FQH effect in groundbreaking work by 
JainEl; the combination of electrons and an even num- 
ber of flux quanta to form these wavefunctions have since 
become known as composite fermions. For odd values of 
m, the wavefunction multiplying the Jastrow factor must 
be symmetric under interchange of two particles ,,-sp that 
such states correspond to composite boson stateS. The 
wavefunctions we study in this work may thus be inter- 
preted as crystals of composite fermions or bosons. 

We therefore propose the following trial wavefunction: 
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Again, 0R i is a single particle wavefunction which is lo- 
calized at R, and lies in the lowest Landau level. Since we 
will consider finite size systems, only those lattice sites 
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within a disk of radius Rd will be included in the set 
of {R;}. In order to make the full wavefunction anti- 
symmetric, we choose to use either a Slater determinant 
(even m) or a symmetric sum of all possible permutations 
(odd to) for the single particle wavefunction part. 

We cannot, however, simply use 4>^f in Eq. (|]^) as 
our single particle wavefunctions, due to the following 
reason. Disregarding the antisymmetrization for the mo- 
ment, the probability density is given by 

i*i 2 ^ni r «-^i 2m ni^^)i 2 ( 2 - 5 ) 

= e 2m E 1 ^ log|r -~ rjl e 2 S l log|0R '' (r,)l . (2.6) 

As in Laughlin's "plasma analogy"!, \^>\ 2 may be thought 
of as the Boltzmann distribution function for a "dual" 
classical system whose effective energy is given by the 
exponents in the above equation, up to an arbitr ary ef- 
fective temperature. The first exponent in Eq. (2.6) is 
identical the 2D (logarithmic) Coulomb energy with cou- 
pling constant m. Each term in the second exponent 
describes an attractive effective potential centered at a 
lattice site R^. Obviously, if the effective Coulomb inter- 
action were absent, the minimum of the effective energy 
would be achieved when = R^ for all i. However, 
due to the effective interaction, the static solution of 
will be moved away from Rj, unless Ri is the center of 
the disk. In general, the electrons will be pushed radi- 
ally away from the center of the disk. Consequently, the 
whole system will spread out and the resulting electron 
density will be smaller than that of the intended lattice. 

In order to prevent this unwanted expansion of the sys- 
tem, for each cj)^ , we will introduce extra zeros ( "ghost 
effective charges" ) outside the physical disk. The ghosts 
are introduced in such a way that if the real electrons 
were fixed at their lattice sites, the total of both real 
and ghost effective charges are symmetrically distributed 
about any given lattice site R^. In other words, the 
ghosts cause each lattice site to look like it is at the 
center of the system by "balancing" out the effective re- 
pulsive force of the surrounding electrons. As a conse- 
quence, each electron will remain centered near its own 
lattice site. Obviously, the specific positions of the ghosts 
depend on the lattice site R^. An example of the way 
ghosts are placed is shown in Fig. H. Although outside 
the physical disk, the ghosts themselves occupy lattice 
sites. Including the ghosts, the single particle wavefunc- 
tion is finally given by 
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where rj^ are the complex coordinates of the ghosts that 
balance out the effective force at R^. 

An interesting property of these wavefunctions is that 
the ghosts may be thought of as "renormalizing" the posi- 
tions of lattice sites. To see this, one can rewrite Eq. (2.7) 
as 
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Once again in the plasma analogy, the logarithms in the 
exponent describe a 2D Coulomb potential caused by ef- 
fective point charges with charge to. Now let us approx- 
imate the point charges {to} by a uniform charge distri- 
bution whose density is the same as the average density 
of the original point charges. To do this, we write the 
sums in the argument of the exponential in the form 

m^2{\og \z-Tjj\+i arg(z - to)} 



to / d 2 -q{\og \z-r)\+i arg(z - rj)}p G (v) (2-9) 



where poiv) is the density of ghost particles. We then 
approximate 

PgM « { P if '^ > R ° ~ ^' < i?G ' (2 10) 

\ otherwise, 

where p is the average electron density, Rd is the ra- 
dius of the physical disk of the finite size system, and 
i?G is the radius of a "ghost disk" , which must satisfy 
i?G > 2i?£>. Provided r is well away from the physical 
disk edge, this approximation should be quite good, and 
we expect corrections to scale as (t/Rd) 2 - Since the real 
part of the integral corresponds to the potential of a uni- 
form charge density toTj in a disk of radius Rq, with a 
circular hole of radius Rd , the real part of the integral 
may be computed using Gauss' law for two-dimensional 
electrostatics. The imaginary part of the integral may be 
computed analytically as well for r <C Rd, yielding the 
approximated wavefunction 



x exp 
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(2.12) 

Note that because the amplitude and phase of <fi' R . have 
been treated on an equal footing, this wavefunction lies 
in the lowest Landau level. Ignoring the unimportant 
constant, </>^. describes an electron in the lowest Landau 
level, centered at a renormalized lattice site (1 — to^)R^. 



Thus the "bare" lattice described by filling 
will be smaller than the real lattice by a factor of 1 - 



R, 



states 
-mv. 
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The physical lattice however is spread back to its original 
size due to the Laughlin-Jastrow correlation. Therefore, 
the above "renormalization" of the lattice compensates 
for the previously mentioned lattice expansion due to the 
Laughlin-Jastrow factor. 

Before we describe the energy calculation for our wave- 
functions, let us briefly discuss about the effect of the 
Laughlin-Jastrow correlation to the characteristics of the 
WC, particularly in connection with the excitation en- 
ergy. In a recent experimental, Jiang et al., have mea- 
sured the temperature dependence of the diagonal re- 
sistance in the reentrant insulating phase slightly above 
v = 1/5. According to their data, the activation gap for 
charge transport is given by E g ~ 0.63K. Surprisingly, 
this energy is much smaller than would be theoretically 
expected. For example, using states of point particles 
whose positions are chosen to optimize the energy, the 
energy to create a point defect such as an interstitial or 
a vacancy has been estimated by many authorsEjIiJ, but 
all the results are an order of magnitude greater than the 
above value of E g . A more recent study of point defects, 
however, shows that the energy can be lowered if the 
Laughlin-Jastrow correlation is introduced between the 
interstitials and the lattice electronsEl Our initial stud- 
ies of interstitial wavefunctions using quantum Monte 
Carlo techniques such as those presented here suggest 
that to reach the very low activation energies seen in 
experiment, one needs to include Laughlin-Jastrow-cpii- 
relations among the ground state electrons as welOE-i 
A discussion of such wavefunctions is deferred to a future 
publication. 



III. COULOMB ENERGY: MONTE CARLO 
SIMULATION 

Since our wavefunction lies completely in the lowest 
Landau level, we only need to minimize the Coulomb 
interaction term in the Hamiltonian. The expectation 
value of the Coulomb energy per electron is written 
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where (• • ■) means the expectation value with respect to 
the wavefunction \& in Eq. (2.4). In the last line, we have 



dropped 1/N and restricted the first integral within a 
single primitive cell at the center of the disk (denoted by 
C) , using the lattice symmetry Since the size of the sim- 
ulated system is inevitably finite, in order to obtain the 
thermodynamic limit, we need to either extrapolate finite 



size results, or use the Ewald sum methodSil. We have 
used the second method in this paper. Details of the cal- 
culation are given in Appendix A, but it must be noted 
here that we have introduced a couple of approximations 
in calculating the Coulomb energy: (I) We have ignored 
the exchange energy, whi ch in practice means that the an- 
tisymmctrization in Eq. (2.4) is dropped. We have tested 
this approximation and find that it is quite good unless 
v is too close to 1/m. The effect of exchange energy will 
be discussed in the next section in more detail. (II) Since 
\^\ 2 when unsymmetrized corresponds to a finite temper- 
ature classical Boltzmann weight, we assume there exists 
a length scale £ c above which fluctuations in the electron 
positions are uncorrelated. We thus use the Monte Carlo 
method to compute the Coulomb interaction between the 
charge density in the central unit cell and the charge out 
to some distance Rs which we presume to be larger than 
£ c . This run is also used to compute the charge den- 
sity in the central primitive unit cell. In order to mini- 
mize boundary effects, we choose Re>, the radius of the 
disk containing all simulated dynamical electrons, to be 
greater than R$, so that electrons close to Rs do not ex- 
perience an environment significantly different than those 
in the bulk. Those electrons between radius Rs and Rjj , 
which are dynamically simulated but not used to com- 
pute the energy, provide an "effective medium" . This 
approach has also-been employed in Monte Carlo studies 
of the FQH effectM 

The interaction of the charge density in the central 
unit cell with charge at distances greater than Rs is com- 
puted by treating the distant charge as static and equal 
to periodic copies of the numerically computed charge 
density in the central primitive cell. This is essentially 
a Hartree approximation. Since this charge density is 
treated as static, one may compute the interaction for an 
infinitely large system using the Ewald sum technique. 
Our method is checked by increasing Rs until the energy 
is unchanged within the errorbars of our Monte Carlo cal- 
culations. Our simulations show that for the wavefunc- 
tion parameters we have studied, £ c is always less than 
4a, where a is the lattice constant. This is also confirmed 
by numerical calculations of individual pair energies, for 
which the result from the simulation is essentially the 
same as the Hartree energy if the pair is separated far- 
ther than 4a. More details of this procedure are discussed 
in Appendix A. 

We have developed a Monte Carlo simulation program 
that computes the Coulomb energy per electron, E c /N, 
using the standard Metropolis algorithmEZl. As a crit- 
ical test of our extrapolation technique, we have used 
our method to compute the energy of the m = state, 
which is identical to the one used in Ref. [52]. Its energy 
can be calculated analytically and our results agree with 
analytic solutions well within the statistical errorbar of 
about 0.05%. The results for more interesting values of 
m are plotted in Fig. |[ Treating m as a variational pa- 
rameter, one can find the value of m which gives the low- 
est energy at a given v. The graph clearly shows that at 
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v = 1/3 and 1/5, the Laughlin state has a lower Coulomb 
energy than any of our wavefunctions. At v = 1/7, how- 
ever, the Laughlin state has a higher Coulomb energy 
than our lowest result. This is consistent with experi- 
ment in that the "true" FQH effect — e.g., vanishing 
diagonal resistivity p xx at zero temperature — has never 
been observed at any inverse odd filling factors below 
v = 1/5. Furthermore, at v = 1/5, the energy of our 
wavefunction is higher than, but very close to that of the 
Laughlin state, which agrees "well with the observation 
of a reentrant insulating phaseEl slightly above v = 1/5. 
This reentrant phase is believed to occur because the pure 
Laughlin wavefunction is the ground state only when v is 
precisely an inverse odd integer. Away from the precise 
filling factors, quasiparticles and quasiholes are present 
in the ground state, increasing the energy. Therefore the 
FQH states have cusps in energy at every inverse odd 
filling factor, allowing the WC state to have lower energy 
in a small but finite range of v right above 1/5. 

Now, let us compare our results with other WC trial 
wavefunctions, particularly the Lam-Girvin forrrO. The 
Lam-Girvin wavefunction also predicts that the phase 
transition from the WC to the FQH effect occurs be- 
tween v = 1/5 and 1/7. As shown in Fig. ^, however, our 
wavefunctions are lower in energy than the Lam-Girvin 
counterpart at all values of v where data is available. 
In other words, our wavefunctions arc closer to the true 
ground state. We believe this difference arises because 
the harmonic approximation neglects rare, but nonethe- 
less important contributions from anharmonic fluctua- 
tions in which two or more electrons come close together. 
In contrast, the Laughlin- Jastrow correlation very ef- 
fectively suppresses density fluctuations at all displace- 
ments of electrons from the lattice sites. This may be 
understood using the plasma analogy for the Laughlin 
statcsu, i.e., the Laughlin- Jastrow correlation is equiva- 
lent to the Boltzmann distribution of a 2D one compo- 
nent plasma (OCP) in which charge density fluctuations 
are suppressed. 

An important difference between the weighting asso- 
ciated with our wavefunction and the Boltzmann weight 
of the OCP is that the electrons are centered at different 
lattice sites in our wave function, while they are centered 
at one single point for the OCP. One of the most signifi- 
cant consequences of this is the following. Let us define 
mo(v) as the value of m for the lowest energy variational 
state at v. Surprisingly, near an inverse odd integer fill- 
ing factor, we find v ~ 1/ (2mo — 1) rather than v *~ 1/ mo 
as in the Laughlin states. For example, at v = 1/7, our 
wavefunction has the lowest energy if m — 4, rather than 
in = 7. Now let us continue to focus on the m = 4 state 
increasing v above 1/7. It continues to be the lowest 
energy state until v reaches about ~ 0.165, where the 
m = 3 state becomes lower in energy. This implies a first 
order transition between the two different m states. This 
phase transition may in principle be detectable in pho- 
toluminescence experiments, although this is presumably 
difficult because the energies of neighboring m states are 



so close together. 

Fig. H shows Ar rms , the root-mean-square value of the 
fluctuation of electrons from their lattice sites. Note that 
Ar rms increases rapidly as v approaches and passes be- 
yond the transition to the m — 1 state. This indicates 
that the single electron probability density becomes less 
and less localized as v increases. However, according 
to the above mentioned plasma analogy, the Laughlin- 
Jastrow factor still tries to force the local density to re- 
main uniform. Therefore, rather than wandering around 
randomly, electrons tend to switch positions and form 
"exchange rings" (Fig. |J) . This means that ring exchange 
energy becomes more and more important. This is most 
easily seen from "snap shots" of the electron configura- 
tion during a Monte Carlo simulation run. We note that 
such ring exchanges are commonly observed in simula- 
tions of melting of the classicaLQCPEJ. In path integral 
descriptions of the FQH effected, coherence among ring 
exchanges play a crucial role in explaining the instabil- 
ity of the WC with respect to a liquid state at v = 1/m 
for small enough m. We believe that quantum coherence 
in ring exchanges may lead to structure in the energy of 
the WC as a function of filling factor even in the insulat- 
ing state, which ultimately could explain the transport 
and photoluminescence anomalies discussed in the intro- 
duction. However, a correct description of this requires 
that exchange be properly included; we therefore defer a 
detailed discussion of this to a future publicationE!]. As 
v approaches and increases past the critical filling fac- 
tor, the exchange rings are observed increasingly often 
in Monte Carlo snap shots. As v increases further and 
Ar rms grows to the same order of magnitude as the lat- 
tice constant, the WC will eventually become unstable, 
giving way to a liquid-like state. This is analogous to the 
melting transition of a conventional solid. In this limit, 
however, the exchange energy is clearly no longer neg- 
ligible and our Monte Carlo analysis ceases to be valid. 
Then, an important question arises: when may exchange 
be ignored? We will address this question in the next 
section. 

Now let us focus on the transitions between different m 
states. First, Fig. || shows characteristics of Ar rms that 
is common to all m. In general, Ar rms is an increasing 
function of v, and as v reaches some point, the system 
undergoes a phase transition to a lower m state. Inter- 
estingly, the values of Ar rms at the critical v is approx- 
imately given by the same value ~ 1.71 regardless of m. 
We believe this implies that derealization plays a crucial 
role determining where the transitions occur. Moreover, 
Fig. H shows that E c /N starts to curve up as v passes be- 
yond the transition value, which is common for all m. We 
believe the derealization and the formation of exchange 
rings are the main reasons for this change of curvature in 
the Coulomb energy. However, it is not yet clear why it 
occurs well below v = 1/m. 
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IV. EXCHANGE EFFECT: PERMUTATION 
MONTE CARLO SIMULATION 

In the previous section, we have seen that the der- 
ealization, which is represented by Ar rms , increases as 
v approaches 1/m from below. Then the exchange en- 
ergy is expected to become more and more important, 
as the overlap of wavefunctions at different lattice sites 
increase. Indeed, we have more or less directly observed 
the formation of exchange rings in the snap shots from 
the Monte Carlo simulations for relatively large values of 
v. Thus, it is clear that as v grows, one must start to 
include the exchange energy in the calculation in order 
to obtain quantitatively reliable results. 

It is very difficult to estimate the exchange energy an- 
alytically for our wavefunction mainly due to the strong 
correlations. However, when m is even, the single par- 
ticle wavefunction part in Eq. (2.4) is a Slater determi- 
nant, and we can take the exchange energy into full ac- 
count by using other Monte Carlo methods such as in 
Ref. ^d]. Our tests with the m = 4 state shows that the 
exchange energy is negligible when v < 1/7. Although 
this method treats the exchange energy exactly, its ap- 
plication is strictly restricted to even values of m, and we 
need to resort to a different method for odd m. 

One way to estimate the relative importance of the ex- 
change effe ct is as follows. The many body wavefunction 
in Eq. (p]j) may be rewritten 



^;) ra EC P fe( r P(i)) 



(4.1) 



where the summation is over all possible permutations 
P. For an odd m, the statistical sign £ p is always +1, 
but for an even m, £ p is either +1 or —1 depending 
on whether P is an even or odd permutation. When 
we ignored the exchange effect in the previous section, 
what we did was to drop all permutations in the above 
summation, except the identity permutation / such that 
= i. In other words, we have approximated the 
above wavefunction with 



^direct = Y[(Zi ~ 2 J ) m 0R I (r l ). 



(4.2) 



Now, we want to define a quantity 7, which measures 
error caused by this approximation, or in other words, 
measures how important the exchange effect is. First, 
the "partition function" maybe written 



Z[S>) = J |*| 2 , 



(4.3) 



where the integral is over all coordinates {ri\. Then, we 
define 
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Z[tf] - Z[V dh 



J2p I Ylijtj 



1 2m 



(4.4) 



(4.5) 



We have used the particle exchange symmetry to reduce 
the number of permutations in each integral from two to 
one. Note that when 7 is small, the exchange effect is 
small and our approximation is good. 

In order to compute 7 numerically, we .have devel- 
oped a "permutation Monte Carlo method"c3, which is 
essentially the same as the usual Monte Carlo simula- 
tion method, except for one important difference: In a 
permutation Monte Carlo simulation, not only the elec- 
tron positions rj, but also the permutation P is treated 
as a configurational variable which is updated, tested, 
and accented (or discarded) according to the Metropolis 
algorithmic!]. Since the integrand in Eq. (4.5) is a com- 
plex quantity, we have separated the phase factor from 
the modulus to sample it. More specifically, we have av- 
eraged the phase factor 



^R,(r t )*^R, (r P(i) ) 
|</>Ri(ri)*<?!>R. ( (rp W )| ' 

using the rest of the factors in the integrand 



m 



(4.6) 



(4.7) 



as the statistical weight in the Monte Carlo simulation. 
We note that this permutation Monte Carlo scheme in 
principle captures the effects of symmctrization of the 
wavefunction exactly. In practice, however, the precision 
of the result depends on how accurately the phase may be 
sampled. As described below, this becomes problematic 
as v — > 1/m. 

The results of the permutation Monte Carlo simulation 
is shown in Fig. [|. It is clear that 7 is negligible if v is 
less than 0.2, 0.16, and 0.14, for m = 3, 4, and 5, respec- 
tively. Beyond those values of v, the exchange effect can 
be significant. The errorbars become quite large as the 
exchange effect gets more and more important. This is 
beca use except for the direct term, the phase factor in 
Eq. (4.6) fluctuates very much, while its average almost 
vanishes. 

Note that we have found earlier from Fig. || that the 
transition between the m = 4 and the m = 3 states oc- 
curs near v = 0.16. Since the exchange effect is negligible 
up to v = 0.16 when m = 4, the m = 4 state is guar- 
anteed to have a lower energy than the m — 3 state if 
v < 0.16, even if the exchange energy is included. Above 
v = 0.16, however, it is currently not known whether 
the exchange effect will raise or lower the total energy 
of the in — 4 state. If it lowers the energy, there is a 
possibility that the transition from m = 4 to m = 3 
actually occur at a higher filling factor than is shown 
in Fig. ||. For m = 5, the exchange effect is negligible 
up to v = 0.14, which is well above v = 0.125 where 
the transition to the m = 4 state occurs. Therefore, for 
this transition, our estimate of the transition filling fac- 
tor is accurate. In principle, however, it is possible that 
a reentrant m — 5 phase occurs within the m = 4 ground 
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state if the exchange effect brings the energy lower than 
that of the m = 4 state above v = 0.14. For m = 3, 
the exchange effect becomes important well before the 
energy level crosses with that of the m = 2 state. How- 
ever, when v > 0.2, not the WC, but the FQH state is 
the ground state, and the energy level crossing between 
different WC states is not physically relevant at T = 0. 
Finally, we note that for m = 3, the exchange terms are 
completely negligible when v < 1/5. Therefore, the com- 
parison of the energy between our correlated WC state, 
the Laughlin state, and the Lam-Girvin state is valid at 
v = 1/5 as well as at 1/7. 



m = 4 states, our estimate of the transition filling factor 
appears to be accurate. 
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APPENDIX A: THERMODYNAMIC LIMIT 



V. SUMMARY 

In this paper, we have studied the correlated WC in a 
strong magnetic field, which is represented by the prod- 
uct wavefunction of the Laughlin-Jastrow factor and the 
Hartree-Fock wavefunction in a triangular lattice. We 
have shown that extra zeros (ghosts) in the single parti- 
cle wavefunction are necessary to balance the expanding 
effect of the Laughlin-Jastrow correlation. 

The energy of the wavefunction has been calculated us- 
ing Monte Carlo simulations and the Ewald sum method. 
Compared to other WC trial wavefunctions in the lowest 
Landau level, particularly Lam-Girvin wavefunction, the 
new state is found to be the lowest in energy. The Laugh- 
lin FQH state has lower energy than our wavefunction at 
v = 1/5 and above, consistent with experiment in that 
the phase transition occurs between v = 1/7 and 1/5. 
The energies of the two states are, however, very close 
to each other near v — 1/5. This is also consistent with 
the experimental observation of the reentrant insulating 
phase slightly above v = 1/5, considering the cusp in the 
FQH state energy due to quasiparticles and quasiholcs. 

Treating m as a variational parameter, we can find 
the value of m which gives our wavefunction the lowest 
energy. Since m takes only discrete integer values, it is 
found that there is a series of first order phase transitions 
between different m states as v changes. For a given 
to, the spatial fluctuations of the electrons grows as v 
increases above ~ 1/(2to — 1), eventually causing the 
energy to curve up. Formation of exchange rings is also 
observed in this limit, which is reminiscent of what occurs 
near melting in classical one component plasmas. 

We have also developed a permutation Monte Carlo 
method in order to estimate when the exchange effect 
becomes important. Our simulation shows that they are 
negligible up to v = 0.2, 0.16, and 0.14 for m = 3, 4, 
and 5, respectively. This ensures validity of our compari- 
son of the energy between different trial wavefunctions at 
v = 1/5 and 1/7. It also shows that when the exchange 
energy is included, the filling factor at which the tran- 
sition between the m = 4 and the m = 3 states occurs 
may be higher, but not lower than our estimate in this 
paper. However, for the transition from the to = 5 to the 



In this appendix, we will explain how one can obtain 
the thermodynamic limit from the finite size simulation 
of the Coulomb e nerg y. We first split the domain of the 
r' integral in Eq. (|3.3[) into two regions 



dr' = I dr' 

s 



dr'. 



(Al) 



The first term is an integral over many, but a finite num- 
ber of primitive cells near the center of the disk, which 
we call the "sampled primitive cells" . In this region, 
the Coulomb energy may be computed directly from the 
Monte Carlo simulations. The second term concerns the 
rest of the whole 2D plane, the "unsampled primitive 
cells" , for which we can not obtain the Coulomb energy 
directly from the simulation. In order to deal with the 
second region, we use the following trick. 

The contribution from the second region may be writ- 
ten 



N 



c 



u 



r — r 



(A2) 



where J c means integral over the central primitive cell 
as in Eq. ( |3.3| ). Now, we assume that the correlation is 
negligible if the separation between r and r' is greater 
than some "correlation length" £ c . If the distance be- 
tween the unsampled primitive cells (U) and the central 
primitive cell (C) is greater than £ c , we may use a Hartree 
approximation and write 



E 



N 



unsmp _ * , , , 

2 Jc J v \r - r> 



(A3) 



where we have defined the expectation value of the local 
density 



p(r) =£<*(r -*<)>■ 



(A4) 



Note that the condition i ^ j in Eq. (A2) is not needed, 



because the domains of r and r' are exclusive so that 
5(r — r, )(5(r' — Tj)Sij = 0. Once the local density profile 
p(r) is obtained in the central primitive cell from the 
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Monte Carlo simulations, the triangular periodicity leads 
to the density profile in the whole 2D plane, and Eq. ( [A3| ) 
can be explicitly computed. 



26 S. Kivelson, D.-H. Lee, and S.-C. Zhang, Phys. Rev. B 46, 
2223 (1992). 



H. Lee, Phys. Rev. Lett. 



Although the integral in Eq. (A3) diverges, this di- 
vergence is unphysical and is easily resolved by recall- 
ing that there is neutralizing background charge in real 
samples. Assuming uniform distribution for the positive 
background charge, the final form of the unsampled part 
is given by 



^unsmp ^ i i 



dr 



,[p(r)-p][p(r')-p] 



(A5) 



where p is the average density of electrons. The above 
integral is now well defined and can be computed using 
the Ewald sum method as in Refs. 33 and Ba. 
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FIG. 1. An example of the underlying triangular lattice 
for a finite size WC. The large circle denotes the bound- 
ary of the physical disk and the filled dots the lattice sites, 
Ri, within the disk. The crosses denote the positions of the 
ghosts, r)j, which balance the effective force at Ri. 
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FIG. 2. The Coulomb energy per electron of the Laugh- 
lin-Jastrow-correlated Wigner crystal as a function of v for 
various values of m. The energy is shown relative to that of 
the uncorrelated Hartree wavefunction (m = 0). The same 
quantity is presented for the Laughlin state (Ref. |l4| ) and the 
Lam-Girvin wavefunction (Ref. Kq) at v — 1/3, 1/5, and 1/7. 
The dashed vertical lines represent the values of v where the 
transition between different m states occur. 
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FIG. 3. Root-mean-square of the electron fluctuation from 
the lattice sites as a function of the filling factor v. The dashed 
vertical lines represent the values of v where the transition 
between different m states occur. 
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FIG. 4. A typical configuration of the three particle ex- 
change ring which can be obtained from a "snap shot" of a 
Monte Carlo simulation. Electron positions are denoted by 
filled circles and lattice sites by empty circles. The arrows in- 
dicate which electron is originated from which Gaussian center 
Ri. 
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FIG. 5. Relative contribution of th e exchange terms in 
the partition function. 7 defined in (4.5) is plotted for several 
values of m. The dashed lines denote errorbars. 
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